Indirect Exchange Interaction between two Quantum Dots in an Aharonov-Bohm Ring 
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We investigate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between two spins lo- 
cated at two quantum dots embedded in an Aharonov-Bohm (AB) ring. In such a system the 
RKKY interaction, which oscillates as a function of the distance between two local spins, is affected 
by the flux. For the case of the ferromagnetic RKKY interaction, we find that the amplitude of 
AB oscillations is enhanced by the Kondo correlations and an additional maximum appears at half 
flux, where the interaction is switched off. For the case of the antiferromagnetic RKKY interaction, 
we find that the phase of AB oscillations is shifted by ir, which is attributed to the formation of a 
singlet state between two spins for the flux value close to integer value of flux. 

PACS numbers: 73.23.-b, 73.23.Hk, 75.20.Hr 



I. INTRODUCTION 

When two magnetic moments are embedded in a metal, 
they induce spin polarization in a conduction electron sea 
and couple each other even they are spatially apart. Such 
indirect exchange interaction, the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction, has been known 
from the 1950s 1]. The indirect exchange interaction 
in magnetic nanostructures is one of basic mechanisms 
for spintronics and it is well understood for ferro- 
magnet and nonmagnetic metal multilayer structures Q . 
However for semiconductor nanostructures, the indirect 
exchange interaction between local spins formed in two 
quantum dots has not yet been observed, in spite of the 
importance as a basic physics and the potential applica- 
tion for semiconductor nanospintronics. 

Recent improvement of the fabrication technique for 
semiconductor nanostructures enables one to make rather 
complicated structures with possibility of the precise con- 
trolling of their parameters. For example, a double quan- 
tum dot (QD) system |4| and the composite system of 
QD and an Aharonov-Bohm (AB) ring have been made 
[a, |(J . The double-dot system was proposed for a can- 
didate of the "qubit" , because in the Coulomb blockade 
(CB) regime a dot with odd numbers of electrons, be- 
haves as a local spin and two dot spins can be entangled 
by introducing the exchange interaction between them 
Q • Such exchange interaction has been also discussed 
from the point of view of the competition between the 
Kondo effect and the antiferromagnetic (AF) interaction 
@ . However the direct exchange interaction was consid- 
ered rather than the RKKY interaction. Investigations 
on the AB ring embedded with QD are aimed at under- 
standing the coherent transport through QD 0, [tj 
and the indirect exchange interaction between two local 
spins has not been addressed. 

So it is intriguing to investigate the RKKY interac- 
tion between two QDs in CB regime embedded in AB 
ring. We will show that the RKKY interaction, the sign 



of which oscillate as a function of the distance (RKKY 
oscillations), is affected by the flux and it dominates the 
transport properties. For ferromagnetic (F) coupling be- 
tween dot spins, the amplitude of AB oscillations is en- 
hanced by Kondo correlations and an additional maxi- 
mum appears at half flux. For AF coupling case, the 
phase of AB oscillations is shifted by tt. 



II. MODEL AND CALCULATIONS 

Figure ^a) shows the schematic picture of an AB ring 
embedded with one QD in each arm. QDs denoted with 
1 and 2 weakly couple to left and right leads. The ef- 
fective Hamiltonian is written with two single-channel 1- 
dimcnsional (ID) leads Hamiltonian Hq and the tunnel- 
ing Hamiltonian Ht as H = Hq + ifp- The lead Hamil- 
tonian is given by H = J2kr=L,Ra £ka\ ra a kra , where 




FIG. 1: (a) Aharonov-Bohm ring embedded with one QD 
in each arm. The system has the parity symmetry along the 
horizontal and vertical axes (dot-dashed lines), (b) The flux 
dependent (left panel) and independent (right panel) particle- 
hole excitation. $ is the flux penetrating the ring. The di- 
rected solid and dashed lines are the particle and hole propa- 
gators, respectively. 
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CfeL(ij)(T is an annihilation operator of an electron with 
quantum number k and spin a in the left (right) lead. 
For simplicity, we adopt the so-called Coqblin-Schrieffer 
model, 



E E 



„T Y n n 
u m(J a 1 a""r'ncr' > 



(1) 



r.r'=L,R' y >< 7 
n=l,2 



=T,I 



for the tunneling Hamiltonian through QDs with odd 
numbers of electrons in CB regime. The Hubbard op- 
erator X™, a = \n, u') (n, a\ describes the spin state of the 
n-th QD and J > is a coupling constant. The an- 
nihilation operator a rna is written using the projection 
(n\kr) of wave function of an electron in the lead r with 
quantum number k at the boundary of the n-th QD, as 

a r ncr=J2k( n \ kr ) a kra EH ■ 

Here we encounter a problem: One needs to know the 
proper wave function including the information on the 
coherent propagation of an electron through arms. Usu- 
ally the scattering theory is suitable for treating the 
electron coherency. However it is complicated to com- 
bine this theory with a theory based on the Hamiltonian 
in the second-quantization representation. In this pa- 
per we circumvent this problem. Rather, we utilize an 
assumption of parity (mirror) symmetry along the hor- 
izontal and the vertical axes (dot-dashed lines in Fig. 
H^a)). Namely, the total Hamiltonian is invariant under 
the interchange of indices L <-> R or 1 2. Such approx- 
imation will make calculations simple and will contain all 
important physics. Any deviations from such a symme- 
try will change the result only quantitatively. 

For practical calculations, it is convenient to introduce 
an annihilation operator of even/odd parity states 0] 
a r ±a = {aria ±a r2 a)/V2, ( a r±cr = J2 k (± I k ) a krcr , where 
|±) = (|l)±|2))/v2 ; we dropped the index r from \kr) 
because of the parity symmetry along the vertical axis). 
Annihilation operators for even and odd parity states are 
orthogonal, 



{a r 



PP'i 



because of the parity symmetry along the horizontal axis. 

When magnetic field is applied, an AB phase factor 
e *0/2 ^ e -j0/2-j mus t kg counted in Eq. QJ, when electron 
tunnels through a QD in the clockwise (anticlockwise) di- 
rection 0[ . The AB phase is written with vector poten- 
tial A as, 



$0 



(2) 



where $o = he/ e is the flux quantum and the line integral 
is performed along the ring in the clockwise direction. An 
AB flux breaks the time-reversal symmetry and it gener- 
ates the main difference between features of the orthodox 
two-impurity Kondo model fll 111 111 ITl ITs| and the 
AB ring embedded with one QD in each arm. Here, we 
note that the magnetic field in leads and QDs is not zero 



for an experiment and it causes Zeeman splitting of elec- 
tron spins. In the following discussions, we consider an 
ideal situation where there is no Zeeman splitting and 
discuss the effect of AB phase on the indirect exchange 
coupling and transport properties. 

The Hamiltonian for the RKKY interaction can be ob- 
tained by the second order perturbation theory in terms 
of J/ef, where ep is the Fermi energy [lLll9|: 



-ffRKKY 



Jrkky(4>) 



The coupling constant Jrkky(</>) can be written as 

J 2 

Jrkky{4>) = y*( 2 + 2cos< ^' 



(3) 



(4) 



where a succeptibility function \ can be found by the 
perturbation theory based on the Keldysh Green function 
technique [l9j |. In the equilibrium, it can be written as 



J 2 X = 7 Re / dede 



j M £ h P ( £ ') - 7p(e)7p( E ') 



4 j e + ir\ — e' 

x {/+(£) -/+(£')} > 



(5) 



where 77 is a positive infinitesimal number and j p (e) — 
JJ2 k (p\k)(k\p)S(e — £fe) is a spectral function of par- 
ity p "electron propagator" . The subscript p represents 
the opposite parity of p, i.e. p = ± for p = Here, 
/ ± (e) = 1/(1 + e ±l3e ) is the electron (hole) Fermi distri- 
bution function, and (3 = T _1 (We use the unit /cb = 1). 
In Eq. Q, a phase dependent factor (2 + 2cos</>) ap- 
pears, which value is related to the four configurations of 
particle-hole excitations - two of which enclose the flux 
[left panel of Fig. IHb)] and pick up a phase factor e 1 ^ or 
e -1 ^ and give term 2cos</>, and the others (right panel) 
are independent of the flux - give term 2. The Eq. (J2J is 
one of the main results of this paper. It shows that by 
means of external flux <f> one can control the amplitude 
of the RKKY interaction but it is impossible to change 
its sign since (2 + 2 cos <f>) > 0. 

Since we consider ID leads, we approximate 7 P (e) as 
for the ID free electron gas with the linearized dispersion 
relation HH: 



l±(e)~J 



1 ± cos 



{M(l + £)} 



(6) 



where fep is the Fermi wave number and I is the length 
of an electron path between two QDs. The argument of 
cosine function is the energy dependent "orbital phase" 
[2pj . i.e. the accumulated phase during electron propaga- 
tion between two QDs. We introduced the cut-off energy 
D = hvpkF, where vf is the Fermi velocity. Here, J 
is written with the density of the states p as J = J p. 
Substituting Eq. © into Eqs. and (J5J, we obtain 
for short distance between two QDs (kpl <§C 2n) the 
ferromagnetic coupling as Jrkky(0) — —2 In 2 J 2 D (2 + 
2 cos 6) and more relevant - the RKKY oscillations of ID 
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free electron gas 1, 21] as a function of / for long distance 
between dots (k F l » 2n) 



Jrkky{4>) 



irJ 2 Dcos{2k F l) 
4U 



(2 + 2COS0). (7) 



Above expressions are obtained by replacing the Fermi 
functions in Eq. JSJ) with those at T = 0, what is valid 
below the characteristic temperature T* defined by, 



_ h hkvi i 

T U Vp 



(8) 



Here, r is the characteristic time scale for an electron 
travels between two QDs. It can be understood from 
the following argument: Electrons deep inside the Fermi 
sea are responsible for the RKKY oscillations. On the 
other hand, electrons with energy e (\e\ <C T*), i.e. elec- 
trons around the Fermi level, are unimportant, because 
such electrons contribute only oscillations whose charac- 
teristic wave length hv-p/T* is much longer than I. Thus 
the RKKY oscillations are insensitive to the temperature 
in the regime T <C T*. However, when the temperate 
reaches T*, the RKKY oscillations are affected by the 
thermal excitations of lead electrons and will be smeared 
out. 

Due to the RKKY interaction, depending on the sign 
of the coupling Jrkky(</>) [Eq. 0], the two dot spins are 
entangled and form a singlet state |0, 0) for AF coupling 
(•/rkky(</>) > 0) or a triplet state |l,m) (m = 0, ±1) for 
F coupling (JrkkyO) < 0). 

In the following, we will discuss how the flux depen- 
dent RKKY interaction modifies transport properties of 
the ring. To calculate the conductance, we adopt the 
third order perturbation theory in terms of J, in or- 
der to take into account the Kondo correlations: They 
are pronounced when the temperature decreases and 
approaches the order of the Kondo temperature Tk ~ 
Dexp(— 1/(2J)). In our case, the RKKY interaction be- 
comes important at the temperature, where the Kondo 
correlations are also unignorable. First we rewrite the 
tunnel Hamiltonian, Eq. ft|. using the vector operator 
of the n-th local spin S n , whose components are defined 



as S'l 



I- 



^U/lT 



and S'2 = (Xft - Xft)/2. Further 

we introduce operators 5 ± = S 1 ± S 2 , which satisfy the 
following commutation relations: 



r cp aPi 



ieijkS k , [Sf,Sj] — ieijkS k 



where is the Levi-Chivita antisymmetric tensor. The 
operator S + does not change the total spin-quantum 
number, while S~ is the operator of the singlet-triplet 
transition. By using abov e op erators, we obtain the sym- 
metrized form of Ht [lj, Ilia . |22| taking into account the 
AB phase as 



p=± I r 
<r,<7'=T,i 

+ cos - (a p RL - S p + v <J Rpa ,5 <r , <r a L p^)+H..c. 
+ isin^ (& RL - S p + v a Rpcr ,6cr'craLpcr' S )+'il-c^ . 



Here, (?„, 



+(-) 



E 



crcr'p r'per 



(10) 



T /0 ! a -/ (T a r p(p) CT denotes effective 



conducting electron spin and is defined with the vector 
Pauli matrix a. Terms proportional to v represent po- 
tential scattering process and for our case, v = 1. The 
first line represents the reflection process and shows that 
the change of parity and the singlet-triplet transition oc- 
cur simultaneously. The second and the third lines de- 
scribe transmission processes. The third line describes 
the singlet-triplet transition without changing the par- 
ity, that is not invariant under the interchange of indices, 
L R or 1 «-» 2, (i.e. the replacement of a T ±a- with 
±o r ±o-). Here it does not mean that the parity symme- 
try is broken: the space inversion transformation changes 
(j) — > — <f>, because it also reverses the direction of the line 
integral in Eq. 

In order to calculate the linear conductance, we adopt 
the diagrammatic technique for the density matrix in the 
real-time domain [23L |2^ | . With the help of the commu- 
tation relations, Eq. the perturbative calculation 
is performed rather systematically (Appendix. El ■ The 
"partial self-energy" which represents the transition rate 
for an electron from the left lead to the right lead accom- 
panied by the triplet-triplet transition Sf^ or preserving 
a singlet state or accompanied by the singlet-triplet 



(9) transition T, L - R (j = 0, 1) is obtained as follows: 
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LR 



^11 



^LR 
J 00 



^LR 



p=± 



3=0,1 




7&(e) lp R {e) Rc <| 1 + £ + ao P ( £ ) + £ 

i=o,i 



T / d£ E « 2 jV £ (e) 7^(e) cos 2 1 + 7 + £ (e) 7p ^) sin 2 | 

/ d£ E 7pi,(£)7pfl(e- A^)Re{l + <Ti p (e) + fTip(e- A 5j )}cos 2 ^ 



(11) 
(12) 



+ 7+ L ( £ ) 7 pfl (e - A j3 -) Re {1 + ^(e) + a p(e - A 3j )} 



sin — 



(13) 



where we neglected the integral Re J ds'j p (e')/(e + if] — 
e'). The subscript j = 1(0) denotes the total-spin quan- 
tum number and j = 0(1). Here A 10 = — Aqi = 
Jrkky(^) and 7^ r (e) = 7 P (e)/ ± (e - /x r ) denotes the 
"lesser" or "greater" Green function where /xl = — /j,r — 
eV/2. The function co(i)p defined by 



co(i)p(e) = de 



e + i?7 — e' 



(14) 



gives the logarithmic divergence related with Kondo cor- 
relations. By substituting Eq. © to Eq. (|14|) . we obtain 



o- ±(e) 
o-i±(e) 



2 J In 



2e 7 P 



ttT ' 
cr ±(e) ± 2 J Re 

Ei (— ifepZ) 



ikvl 



in 



2T* 



for V = and e < T,T*. Here Ei(x) denotes the ex- 
ponential integral function and 7 s» 0.577 is the Eulcr 
constant. Equation 10} supplemented with Eq. J5j and 
Eqs. I|ll|) . i|12|l and (|13|l are main results of this paper. 

Using the partial self-energy, Eqs. fTTft . (JT^J and (fT5|>. 
the current can be expressed as 



i,i'=o,i 



P, {S^-(L^i?)}, (15) 



where probabilities Po for a singlet state and Pi for each 
of particular triplet states (we consider no Zeeman split- 
ting), can be obtained for the linear response from the 
Boltzmann distribution as 



Po = 



1 



Pi = 



1-Po 



(16) 



1 + 3 cxp(-/3J RK KY (</>)) ' 3 
The linear conductance is defined as G = limy^o dl/dV 



III. RESULTS AND DISCUSSION 



For the AB ring geometry without quantum dots in 
arm, the conductance oscillates as a function of the flux <p 
fl2| . Furthermore, because of the orbital phase, the con- 
ductance also oscillates as a function of the length of the 
arm I for enough low temperatures: As the thermal exci- 
tation of lead electrons scrambles various orbital phases, 
the oscillatory component would be reduced for temper- 
atures above the characteristic temperature T*, where 
characteristic length hvp/T reaches /. When one local 
spin, i.e. QD, is embedded in each arm of the AB ring, in 
addition to the oscillatory component, the non-oscillatory 
background of oscillations related to spin-flip processes 
will appear: Spin-flip processes do not contribute to the 
interference effect [25j because if a local spin is flipped, 
we can determine the path which an electron propagated. 
Such non-oscillatory background reduces the portion of 
oscillatory component. However, if we take account of 
the RKKY interaction, the oscillatory component can be 
enhanced from the following mechanism: First, according 
to Eq. (|16|l . the probabilities for singlet Po and triplet Pi 
states will be affected by the RKKY coupling constant 
Jrkky (</>), which is a oscillatory function in terms of 4> 
and I, Second, as the conductance would be sensitive to 
a state of local spins, it would show also the oscillatory 
behavior related to the oscillations of Jrkky (</>)• Such 
RKKY dominant oscillations one could expect for the 
enough low temperature T<C |./rkky(0)| <CT*. 

In the following, we will discuss the properties of our 
system for temperatures where the thermal scrambling 
of orbital phases is unimportant T <C T* and above 
the Kondo temperature T ^> Tr. We note that as 
I Jrkky (0)1 ^lT*, the modification of orbital phases by 
inelastic spin-flip scattering events is also unimportant. 
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A. / - dependence 

First we will discuss the RKKY oscillations without 
magnetic flux <j> = as a function of the distance I be- 
tween two QDs. Figure Eta) shows the RKKY oscilla- 
tions of the coupling constant Jrkky(O) as a function 
of the length of an electron path I. It oscillates with 
the period of hpl/it — 1, and shows local minima at 
integer values of k-pl/Tr corresponding to ferromagnetic 
(F) coupling and local maxima at half-integer values of 
k-pl/n corresponding to the antiferromagnetic (AF) cou- 
pling between the spins. The amplitude of the oscillations 
decay with 1/kpl as predicted for RKKY interaction in 
quasi ID geometry |2l[. In Figure E^b), there is a plot 
of the probability Pq of the singlet state. At low tem- 
perature T < | Jrkky(0)| (the solid line), a singlet state 
(triplet state) is formed when value of kpl/ir is close to 
half-integer (integer). As the temperature increases (the 
dashed line for T ~ | Jrkky(0)| and the dotted line for 
T ^S> | Jrkky(0)|) the amplitude of oscillations is sup- 
pressed and system approaches uniform distribution be- 
tween the singlet and triplet states - Pq = Pq = 1/4. 
There are also the oscillations of the conductance [Fig. 
Etc)] with the period of fepi/vr = 1. In the same way 
as in Fig. E^b), the amplitude of oscillations is sup- 
pressed for T ^> |Jrkky(0)|, what indicates that in 
regime T < | Jrkky(0)| the conductance oscillations are 
mainly determined by the RKKY interaction. Experi- 
mentally, it can be difficult to control the length of arms 
keeping other parameters fixed. However, the conduc- 
tance oscillations would be possible to observe by chang- 
ing the Fermi wave number ftp, by controlling the carrier 
density of 2D electron gas (2DEG) with an additional 
gate. 



2x10" 



B. 



dependence 



Though above discussions suggest that the RKKY in- 
teraction dominates the length dependent conductance, 
it would be more convenient experimentally to measure 
the flux dependence. In the following wc will discuss the 
modification of AB conductance oscillations by the pres- 
ence of RKKY interaction. 

As we mentioned before by means of external flux cj) 
one can change the amplitude of the RKKY interaction 
but not its sign since (2 + 2 cos <ft) > 0. In the particular 
experimental situation depending on the length I of the 
arm and Fermi wave vector fcp the spins can be coupled 
ferromagnetically or antiferromagnetically. By means of 
flux <j) one can control the strength of the interaction 
but does not switch between them. For this reason it is 
generic to discuss three typical situations, for which we 
are able to get analytic results. These three cases are 
classified by the value of the RKKY coupling constant: 

(i) the uncorrelated local spins case (| Jrkky(</>)| <?"), 

(ii) the ferromagnetic coupling case (— Jrkky(0) 3> T) 
and (hi) the antiferromagnetic coupling ( Jrkky(</>) ~^>T). 
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FIG. 2: Length dependent (a) RKKY coupling constant, (b) 
probability for singlet state and (c) conductance for T/D = 
5 x 10" 5 (solid line), 10~ 4 (dashed line) and 10~ 3 (dotted 
line). Parameters are taken as <j) — and J = 0.04. 



(i) Uncorrelated local spins limit is realize for high tem- 
perature | Jrkky(^)| <C T or for the flux w ir + 2mi 
since then, according to Eq. 10}, the RKKY interaction 
is weak, Jrkky(0) - * 0. In this case (i), the local-spin 
state is distributed with equal probability among a sin- 
glet state and a triplet state P\ » Pq rs 1/4 [see Eq. l|To|l]. 
The conductance is expressed as 

G 



— ~ {kJ) 2 \ v 2 (l + cos0 cos 2 fc F ;) 



3 1 + 4Jln 



2e 7 £>\ 
ttT J 



(17) 



where Gk = e 2 /h is the quantum conductance. The first 
term, which is proportional to v 2 and thus independent 
of spin-flip processes, is attributed to the phase coher- 
ent component of the cotunneling process. It shows the 
ordinary AB oscillations. The second term in Eq. I|17jl. 
which is related to spin-flip processes (does not depend 
on (ft), form the background of AB oscillations. We can 
see that with decreasing of temperature the Kondo cor- 
relations enhance the background: The second term can 
be interpreted as the parallel conductance through two 
independent spin- 1/2 local moments whose conductance 
is enhanced by Kondo correlations [2(j. In the third or- 
der contribution in J in Eq. I|17[l . there is no interference 
related to the orbital phase k-pl, which was pointed out 
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by Beal-Monod 0. We explicitly showed by Eq. i(TT|) 
that there is also no interference related to the AB phase 
in the third order contribution. 

(ii) The ferromagnetic coupling for — Jrkky(</>) 3> T : 
In this case, two local spins form a triplet state P\ m 1/3 
and Pq » [see Eq. Thus, the conductance is that 

of S = 1 Kondo model plus the potential scattering. For 
the case of long distance between QDs (kpl 3> 1), 



£ * 2( 7 rJ) 2 



- 2 o 2T* 

4 J cos fcp^ cos — In — — 
2 7rT 



v 2 - 2e 7 Z?\ 
1 + — + 2 J In — — J (l + cos (/> cos 2 fcpZ) 

(18) 

For the opposite case, kpl 1, we obtain the same equa- 
tion as Eq. I|18[) with replacing T* in the logarithm by 
e~~<D. The striking feature is that as opposed to the 
case (i), the Kondo correlations enhance the oscillatory 
component as it is shown in the second term of Eq. (|18|l . 
Loosely speaking, two spins are no longer independent 
phase-breaking scatterers because they "observe" each 
other and the Kondo correlations enhances the AF cou- 
pling of each QD spin to the conducting electrons spins. 
The first term of Eq. i|18|l shows the logarithmic diver- 
gence, whose the cutoff energy is equal to the character- 
istic temperature of the orbital phase coherence T* . This 
term appears because the spin-1 moment stretches over 
Using Eq. I|18[) . we can relate the F coupling of spins by 
RKKY interaction [Fig. with the maximum in the 

conductance [Fig. Hfc)] around integer values of kpl /it. 

(hi) The antiferromagnetic coupling, Jrkky(</>) T: 
In this case two local moments form a singlet state Pi « 
and Pq « 1 [see Eq. Ijl6jl ]. As the singlet state is decou- 
pled form lead electrons, i.e. electrons flowing through 
QDs cannot excite local spins to a triplet state, so only 
the potential scattering process contributes to the con- 
ductance: 



G/G K ^ (ttJ) 2 w 2 (1 +cos0 cos 2 k F l). 



(19) 



Because we consider the Coulomb blockade regime, the 
cotunneling current is very small. It is the reason why the 
conductance is suppressed around each half integer value 
of kpl /tt [Fig. EIc)], where the RKKY coupling is antifer- 
romagnetic [Fig. HI a)]. Here we note that situation (ii) 
and (iii) are not realized for the flux (j> w tt + 27rn since 
•/rkky(^) is small there and the limit (i) is approached. 
For | Jrkky(0)| > T by means of the flux we can tune 
between (ii) and (i), or (iii) and (i) but not between (ii) 
and (iii) situations. 

For above three cases, we obtained simple analytic re- 
sults and clarify that the local-spin state due to RKKY 
interaction causes the pronounced effect on the conduc- 
tance. Next, we will analyze the conductance of the sys- 
tem for the full range of the flux <j> and discuss the ad- 
ditional structures caused by the flux dependent RKKY 
interaction, which can be an evidence of the RKKY in- 
teraction in our system. Figures |3Jl-a) and (2-a) show 
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FIG. 3: Flux dependent RKKY coupling constant ((1-a) and 
(2-a)), probability for singlet state((l-b) and (2-b)) and con- 
ductance ((1-c) and (2-c)) for J = 0.04. Panels (1-a), (1-b) 
and (1-c) correspond to the F coupling case (hpl/fin) = 5) 
and panels (2-a), (2-b) and (2-c) correspond to AF cou- 
pling case (kpl/(2n) = 5.25). Flux dependent conductance 
for (3) M/(2vr) = 50 and (4) 50.25. The solid, dashed, 
dotted, and dot-dashed lines show the results for T/D = 
5 x 10~ 5 , 10~ 4 ,2 x 10~ 4 , and 10~ 3 , respectively. 



the RKKY coupling constant Jrkky(</0 as a function of 
the flux (f>/(2ir). The former shows plot for F coupling 
case (kpl/ir is an integer) and the latter show the plot for 
AF coupling case (kpl /tt is a half-integer). The panels 
(1-b) and (2-b) are the corresponding plots of the prob- 
ability for the singlet state for various temperatures and 
the panels (1-c) and (2-c) are plots of the conductance. 
In the vicinity of zero flux 7r = 0, electron wave functions 
constructively interfere and thus the maximum RKKY 
interaction is induced [panels (1-a) and (2-a)]. For F 
coupling case, a triplet state is formed, i.e. Pq ~ 0, at 
low temperature [panel (1-b)] and thus the conductance 
is enhanced [panel (1-c)] as discussed in (ii). For AF cou- 
pling case at low temperature, a singlet is formed [panel 
(2-b)] and the conductance is suppressed [panel (2-c)] as 
it was discussed in (iii) . At half flux, electron wave func- 
tions destructively interfere and the RKKY interaction is 
switched off [panels (1-a) and (2-a)]. Surprisingly at half 
flux we can observed the maximum in the conductance for 
both situations F and AF. According to discussion in (i), 
this maximum is caused by the term in Eq. (|17|) . which 
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does not depend on the flux, and which corresponds to in- 
coherent transport thought the two independent spin-1/2 
local moments related to Kondo correlations. Especially 
for AF coupling case, it leads to the effective phase shift 
of AB conductance oscillations by ir [panels (2-c)]. 

In order to compare our results with the limit, where 
the RKKY interaction is negligible, we show curves of 
AB oscillations for | Jrkky(</>)| C T in panels (3) and 
(4) . As discussed in (i) , the component of the ordinary 
AB oscillations is very small. The Kondo correlations 
only enhance the background and they do not promote 
characteristic structures as the case of the AF or F cou- 
pling. 

Here we note some features to distinguish experimen- 
tally the RKKY dominant oscillations from the ordi- 
nary AB oscillations. The first feature is the charac- 
teristic temperature below which the oscillations can be 
observed: The characteristic temperature of the ordi- 
nary AB oscillations T* is higher than that of RKKY 
dominant oscillations |Jrkky(0)| by the factor ~ J -2 . 
One can point out that the RKKY dominant oscilla- 
tions is sensitive to the temperature. The second fea- 
ture is the temperature dependence of the amplitude of 
oscillations: Suppose we decrease the temperature from 
enough high temperature T 3> |</rkky(0)|, where sin- 
glet and triplet probabilities are Pj w Pi « 1/4 and 
the conductance is expressed by Eq. 1|17|) . As tem- 
perature is lowered, singlet and triplet probabilities are 
modified as P ~ 1/4(1 - 3JrkkyO)/(4T)) and Pi » 
1/4 (1 + Jrkky(^)/(4T)). Therefore, the correction de- 
pending on both the orbital phase and the AB phase 

f / - 2e 7 D \ 
SG ^-(,jf 3 l+4Jln— (2 + cos 2 Mcos0) 

Jrkky(^) 
4T ' 



2u 2 (l + cos 2 Mcos( 



emerges. It grows as oc (In T)/T\ the logarithmic correc- 
tion is related to the Kondo correlations. We expect that 
with the help of the Kondo correlations one can distin- 
guish the RKKY dominant AB conductance oscillations 
from the ordinary AB oscillations. 

Here we will note on the Onsager symmetry. For the 
two-terminal geometry, it means that the conductance 
is an even function of the flux. We can see that the 
RKKY coupling constant Jrkky(</>) is an even function 
of the flux (Eq. (0J). This property depends only on 
the symmetry of the Hamiltonian under the inversion of 
time and magnetic field an d does not depends on the 
assumption of the mirror symmetry. 

Finally we discuss on the range of parameters. For a 2 
DEG system at an AlAs/GaAs heterostructure, the car- 
rier density of which is typically w 3.8 x lO 15 m _2 0, the 
Fermi energy and the Fermi wave length are £f ~ D rs 
14meV and 2n/k F w 40 nm, respectively. The RKKY 
coupling constant Jrkky(0)~J 2 P?/ (&f0 should be larger 
than the Kondo temperature Ir ~ Z)exp(— 1/(2 J)), 
I Jrkky(0)| 3>Tk, otherwise, each spin-1/2 local moment 



forms Kondo singlet and are screened out and thus the 
RKKY interaction is unimportant. In our calculations, 
we put J = 0.04 which gives the small Kondo temper- 
ature, T K ~ 3.7 x 10~ 6 P> < |Jrkky(0)|. Because we 
adopted the perturbation theory, the temperature should 
be above the Kondo temperature, T 3> Tk- In order to 
obtain the large RKKY interaction Jrkky(O) > T we 
put the size of the ring as kpl w 5 x 2ir (about 200 nm). 

For the case of a small AB ring, the Zeeman split- 
ting could become important. In order to reduce the 
Zeeman splitting Ez keeping the number of fluxes con- 
stant, one could increase the size of the AB ring be- 
cause Jrkky(0)/-Ez ~ I- In such a case, the consider- 
ation of the Kondo regime T < Tk, could be needed 
because RKKY interaction was also reduced so the limit 
Tk ^ I <Prkky(0)| was approached. 



IV. SUMMARY 

In conclusion, we have theoretically investigated the 
RKKY interaction acting between local spins, i.e. two 
QDs with odd numbers of electrons in CB regime, embed- 
ded in the AB ring. We assumed the parity symmetry of 
the system and such an assumption does not change the 
result qualitatively. We calculated the RKKY coupling 
constant and the conductance above the Kondo temper- 
ature, T > Tk. The RKKY coupling constant, the sign 
of which oscillates as a function of the distance, also de- 
pends on the flux and the distance between two QDs. 
When the RKKY interaction is ferromagnetic, two local 
spins form a triplet state around zero flux, where the 
electron wave constructively interfere, and thus the max- 
imum RKKY interaction is induced. As the temperature 
decreases, the amplitude of AB oscillations is enhanced 
by Kondo correlations, which is the distinctive difference 
between the ordinary AB oscillations and those of the 
ferromagnetically coupled two local spins. The maxi- 
mum was found at half flux where the RKKY interaction 
is switched off and the conductance is described by the 
parallel conductance of two independent spin-1/2 local 
moments whose conductance is enhanced by Kondo cor- 
relations. When the RKKY interaction is AF, the phase 
of AB oscillations is shifted by 7T. It is because around 
zero flux, where we obtain the maximum AF interaction, 
two local spins form a singlet state, which is decoupled 
from the lead electrons. 
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APPENDIX A: THIRD ORDER PERTURBATION 
THEORY 

In this Appendix, we present detailed calculations of 
the third-order partial self-energy in terms of J on the 
basis of the diagrammatic technique in the real-time do- 
main [23l,l24l |. Figure 0] shows the second order diagrams 
for the partial self-energy representing the transition pre- 
serving the total-spin quantum number j, Ej^ ((a-1), 
(a-2) and (a-3)) and the singlet-triplet transition, Ejj^ 
(j = 0,1) ((b-1) and (b-2)). Green functions of lead 
electrons are represented by directed solid lines, which 
are also called "reservoir lines", and solid lines on the 
Keldysh contour (two horizontal lines) represents propa- 
gators of local spins. Here diagrams (a-1) and (a-3) rep- 
resent different processes. For the former case, we must 
count factor —1 for the vertex denoted with Sj when 
(j = L We omitted diagrams which could be obtained by 
applying the mirror rule psf . 

Following the rules in Ref. [2^, the diagram (a-1) plus 
its mirror diagram can be calculated as 



(a-1) 



(a-2) 



(a-3) 



J.m 



j,m + \ 



J.m 




(b-1) 



(b-2) 



J.m 



j,m + \ 




FIG. 4: The diagrams for the second order partial self- 
energy representing the transition preserving total-spin quan- 
tum number j, [(a-1), (a-2) and (a-3)] and the singlet-triplet 
transition [(b-1) and (b-2)]. Directed lines represent propa- 
gators for lead electrons. Thick solid lines on the Keldysh 
contour (two horizontal lines) represent propagators for the 
local spins. 



utilizing the commutation relations Eq. @ as 



(a-1) 



P =± 
-j<m<j 



V / d£ 7+L( £ b pfl (e)cos 2 - 



x Re (j, m| 25+ \j,m) 



(Al) 



The results for the diagrams (a-2) and (a-3), which 

we term E£- a ^ and S^ ( ' a ^ , can be obtained from 

Eq. (SB by changing 2S+ to S"+S'+ and to 2v 2 , respec- 
tively. In the same way, the diagram (b-1) plus its mirror 
diagram is calculated as 



(b-i) 



E 



P =± 

-j<m<j 



de 7+ L ( e )7p fl (e- A,j)cos 2 7 



RcO',m|2S- \j, 



(A2) 



where A^-j is the energy difference between the total-spin 
quantum number j state and j state. The result for the 
diagram (b-2), which we denote by ^ff' b 2 , is obtained 
from Eq. JUJ by replacing of 2S~ 2 by 5+5+. 

The third order diagrams give the vertex correction to 
the second order diagrams. Figures Ufa'-l), (a'-2), (b'-l) 
and (b'-2) show the correction for the vertex on the upper 
branch of diagrams (a-1), (a-2), (b-1) and (b-2) in Fig.0J 
respectively. Except for the topmost two diagrams, we 
omitted the lower branch of each diagram, which is ex- 
actly the same as for the corresponding diagram in Fig. 
QJ The left diagrams and the right diagrams show direct 
tunneling processes and exchange processes, respectively. 
We did not show the correction for the diagram (a-3) be- 
cause it is proportional to j< m <,-(j, w|S'+|j, m) and 
thus vanishes. For example, the topmost two diagrams 
in Fig. [3] plus their mirror diagrams are calculated by 



^Lft(si-l) correction 

2i 
43 



Im / deide2de3—^- 

rp ^ 



lLp( £ lhRp( £ 3) 2 4> 

COS — 

£3 + 17] 2 



+ 



I £1 - £ 2 + ir\ £2 - S3 + it] J 

,2 



J de 7^(e)7 fl ,p(e)cos 2 



n TU <^4^am\2Sf\j im )\, 



(A3) 



where o\ v is defined in Eq. I|14|) . Here we counted the 
minus sign for a loop with three vertices in the anticlock- 
wise direction and we dropped terms except for the renor- 
malization of the transmission probability. We checked 
that terms which we dropped are canceled out by the 
other diagrams than those depicted in Fig. |5j i.e. dia- 
grams in which the position of a lower vertex is inbetween 
upper two vertices. Further we neglected the integral 
Re J de'j p (e')/(e + ir\ — e') which is at most ^Je/D for 
k F l<s:l or ~J{\e\/D + l/(k F l)} for fc F Z>l. By adding Eq. 
(|A1I) and Eq. (|A3|> . we obtain Eq. ljAl|) with replacing 
2S+ 2 by2(l + a lp (£)/4)S+ 2 . 

Other two diagrams in panel (a'-l) can be calculated in 
the same way. In Fig. [5] we omitted diagrams obtained 
by reversing direction and indices for spin and lead of 
reservoir lines. By calculation of all such diagrams and 
adding them to Eq. (|A1|I , we obtain Eq. (|A1|) with replac- 
ing 2S^ by 2z'j J Sf , where the renormalization factor 
is given by 



"jj 



= 1 + \ Lp(e) + ^-Aa)W e -A„)- 
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FIG. 5: The third order diagrams: Each four diagrams of 
(a'-l), (a'-2), (b'-l) and (b'-2) show corrections for the vertex 
on the upper branch of the diagram (a-1), (a-2), (b-1) and 
(b-2) in Fig. 0] respectively. 

In Fig. El we did not show the lower vertex corrections, 
which are given in the same way as the upper vertex 
corrections. By counting lower vertex corrections, z'-- is 
modified as 

^ = i + ^ (e) + ^M+^zM , ( A4) 

Finally the third order contributions change Ej^ a to 
X;^ R ^ a ~ 1 , where the latter is obtained from the for- 



mer by replacing 2S+ with 2zjjSf . For the other 
diagrams than those of panel (a'-l), we can repeat the 
same discussions as above. The result for diagrams (a-2) 
and (a'-2) and their derivative diagrams, which we term 

£|f (a " 2)(2) , is obtained from Eq. (HJ by replacing 2S+ 2 
with ZjjS±S^. 

By calculating diagrams in (b'-l) and their derivative 
diagrams, and adding them to the diagram (b-1), we ob- 
tain E 1 ^ 13 , which is the same expression as Eq. (IA2|) 

2 2 

with replacing 2S~ by 2ZjiS~ . Here 

Zfj = 1 + ai p (e) + <r lp (e - A- }j ). (A5) 
The result for the diagrams (b-2) and (b'-2) and their 

derivative diagrams, which we term Y> L -P^ h 2 ^ 2 ^ is ob- 

tained from Eq. (|A2f) by replacing 2S~ with ZjjS±S~. 

tt ii i • • •nrR( a - 1 )( 2 ) v^f»(a-2)(2) . 

Finally, by summarizing Eff y A Y.j } m ' and 

T,ff- ^ , we obtain the first term of Eq. fTTJl for j = 1 
and the first term of Eq. i|12|) for j = 0. By adding 

S £R(b-l)(2) E i iJ (b-2)(2) we fi tcrm of 

For now, we have explained only the diagrams related 
to the time-reversal symmetric term, corresponding to 
the first and second lines in Eq. I|10|l . Diagrams related 
to the time-reversal symmetry breaking term, the third 
line in Eq. I|10(l . are obtained from Figs. 0] and [5] by 
changing the parity indices of the right reservoir lines. 
For example, the corresponding diagram of (a-1) is cal- 
culated as 

p=± 

m=±l,0 

x Re(j,m\2Sf\j,m). 

For vertex corrections, the change in the parity indices 
of the right reservoir lines corresponds to the operation 
of the replacement of <j\ v by ao p . Thus the second terms 
of Eqs. i|ll[l . (|12fl and (|13|l can be obtained from first 
terms by replacing <Ji p , cos 2 (</>/2) and 7^)^ with crop, 
sin 2 ((/>/2) and l p{p)w respectively 
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